Multisymplectic Geometry and Multisymplectic Preissman 
Scheme for the KP Equation * 



fsj i Tingting Liu, Menzhao Qin 

O ; CCAST (World Laboratory), 

^ ■ PoBox 8730, Beijing 100080, China 

^ . Institute of Computational Mathematics, 

d I Academy of Mathematics and Systems Sciences, 

Chinese Academy of Sciences, 
PoBox 2719, Beijing, 100080, China 
ttliuOlsec .cc.ac.cn, qmzQlsec .cc.ac.cn 

> 

o 

o . 

' Abstract 
O. 

CN| , The multisymplectic structure of the KP equation is obtained directly from the 

' variational principal. Using the covariant De Bonder- Weyl Hamilton function the- 

, ories, we reformulate the KP equation to the multisymplectic form which proposed 

CIh' by Bridges. From the multisymplectic equation, we can derive a multisymplectic 

numerical scheme of the KP equation which can be simplified to multisymplectic 
forty-five points scheme. 



> 
X 



1 Introduction 



^ ' The generalized Kadomtsev-Petviashvili (GKP) equation is 

(ut + dxfiu) + Uxxx)x + cruyy = (t > 0, -oo < x,y < oo) (1.1) 

where cr is a constant, f{u) is some smooth function. The usually KP equation occurs for 
f{u) quadratic and it is regarded as a two-dimensional generalization of the Korteweg-de 
Vries (KdV) equation. It describes the evolution of long water waves of small amplitude 
if they are weakly two-dimensional. In the case of f{u) = Su^ and cr = — 3, Equation 
(1.1) is usually called the KPI equation, whereas the KPII equation with f{u) = 3n^ 
and (7 = 3. As a soliton equation important from analytical and numerical point of view, 
the KP equation is one of the few known completely integrable equations in the multi- 
dimensional soliton equations. Thus, in the last few years, considerable attention has 
been paid to the KP equation. Although considerable interest has been focused on the 
KP equation, the numerical scheme analysis literature for the KP equation is extremely 
small. As far as we are aware of, Katsis proposed explicit finite difference method [15], 
the results of evolution of lump solution for the KP equation was given by A. A. Minzoni 
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[17], X.P. Wang ct. studied the instability of a generalized KP equation [16], B.F. Feng 
and T. Mitsui took the linearized implicit method to the KP equation [13]. 

In this paper, we try describe KP equation in the language of multisymplectic geom- 
etry. Recently, for first order field theory, i.e., the Lagrangian density depends on the 
state variables and their first order derivatives, Marsden, Patrix and Shkoller [7] derived 
numerical methods for the first order field theories. However the lagrangian density of 
the KP equation is not first-order, therefore MPS theory can not be applied directly. In 
[10], the authors focus their attention on the KdV equation whose Langrangian density 
is second order. The langrangian density of the KP equation is truly third-order. In this 
paper, we give the multisymplectic structure of the KP equation directly from the vari- 
ational principal. In [5], the author proposed the Cartan form is not necessarily unique, 
we find it was caused by higher-order mixed multiple integral in using Stokes' formula in 
actually calculus variation. 

In Lagrange mechanics, we know the Euler-Lagrange equation can be write as 

(1.2) as regular Hamilton equation 



Taking the Legandre transform of Lagrange density L : = |^ , we can rewrite equation 



dp'- _ dH 
dt dq" 

dq' ^ dH 

dt 9»» 



(1.3) 




where H = p^q^ — L{q^,cf,t). With the covariant De Bonder- Weyl Hamilton function 
theories [18], we can reformulate the partial differential equation to the following form 

where 7rf = gf^- According to this methods, we can rewrite KP equation to the 
multisymplectic form that introduced by Bridges [3] . 

Multisymplectic equation have the important multisymplectic conservation laws. In 
the numerical study, we also hope that the numerical approximations can preserve the 
multisymplectic conservation laws. Similar to the method [9], we show that the Preiss- 
man scheme is a multisymplectic scheme for the KP equation. Though the Preissman 
scheme is multisymplectic, it need more computational memory, so we reduce it to a mul- 
tisymplectic forty-five points scheme. Using the forty-five points scheme, we get some 
numerical results on soliton and solitary waves over long time intervals. 

In section 2, we describe the multisymplectic geometry of the KP equation entirely 
in the framework of the variational principal. Section 3 is devoted to the analysis of 
multisymplectic Preissman scheme and reduce it to a multisymplectic forty-five points 
scheme . In section 4, some numerical results on soliton and solitary waves over long 
time intervals be given. 
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2 Multisymplectic Geometry of the KP Equation 



We now review some aspects of multisymplectic geometry. 

Let X be an orientable (n + l)-dimensional parameter space ( which is usually space 
time ) and let ttxy ■ Y ^ X he a fiber over X. Section (p : X ^ Y of this covariant 
configuration bundle is the physical fields. Coordinates on X are denoted by x^,fj, = 
1, 2, • • • , n, 0. In general, denotes the time coordinate. The parameter n denotes the 
number of spatial variables. In this paper, we just discuss the case of n > 0. Adapted 
coordinates on Y are along the fibers 1^ = tt^y{x),x G X , A = 1,2, ■ ■ ■ , N . N 
denote the fibers dimensions. Consider a k^^ order lagrangian density L, viewed as a 
fiber-preserving map L : J^Y — > A"'"'"''^^. J^(Y) denote the A;*'*-order jet bundle over Y 
which can be induced by J^(- • • (J^(y))). We let T^X denote the tangent space of X at 
X. Denote the derivative of the map ttxy in the direction w by T^^xy ' 

At first, we introduce the first jet bundle. 
Definition 2.1 The first jet bundle over y is a fiber bundle denoted by J^{Y) whose 
fiber over y E Y^ = 7r^y(x),x G X consists of those linear mappings 7 : T^X —>■ TyY 
satisfying 

Tt^xy 07 = Identity. 

If 99 : X — > y is a section of Y, j^{ip) is a section of J^{Y) and in coordinates, j^i^p) is 
given by 

{x^^, ip^{x^'),d^ip^{x^)), M = 1, 2, • • • , n, 0, 

Similarly, higher order jet bundle J"^{Y) is defined by J^(J™^^(y)). 
Definition 2.2 The A;*'*-order jet bundle over y is a fiber bundle denoted by J^(y) whose 
fiber over 7 G j''~^{Y)y,y G Y consists of those linear mappings s : T^X — Tjj''~^{Y) 
satisfying 

Tjr^ /fc-iy ° s = Identity. 

We let j^{ip) = ■ ■ (j^(<^)) : x Txj^~^{'p) denote k^^-oxdex jet prolongation of the 
section : X ^Y,'m. which j^~^{if) is a section of jet bundle J^^^i^). Thus, is 
given in coordinates 

(x^ ^\x^), d,^\x^^), • • • , d,^d,^_, ■ ■ ■ d,,^^{x^)), 

Given a k^^-ordei lagrangian density L : J^Y A"'*"^,t, the basis geometric object in 
the classical calculus of variations is the (n + l)-form 6l on J^'^~^(y), which was called 
the Cartan form. 

The KP equation can be written as 

{2ut + 6uux + Uxxx)x + cruyy = 0. (2.1) 

The 2 multiplying ut is added for notational convenience, it can be eliminated by scaling 
t. In this paper, we consider the KPI equation. 

To put the KP equation in the variational frame work, we let ipxx = u, then (p satisfies 
equation: 

"^Vxxxt + ^fxxVxxxx + 

XXX "I" Vxxxxxx "I" (^^xxyy 0- (2-2) 
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The search for a variational principle is equivalent to the inverse problem of the calculus 
of variations, i.e., the existence and formulation of variational principles for systems 
of nonlinear partial differential equations. The existence of a variational principle for a 
differential equation is equivalent to determining whether or not an operator is a potential 
operator. According to Vainberg theorem [2], in order that operator A'^ be potential 
operator, to summarize, it is necessary and sufficient that the Gateau derivative of the 
operator A*" is symmetry. The theorem stated in terms of the Gateau derivative, but we 
assume that the Prechet derivative exists in application. Let N is an operator which 
defined in an appropriate function space E (typically E is Banach space), then N be 
potential operator if = N'^, where 

N^ip = hm ^ ^ ^ = [— Ar(« + e^)] |,=o, 

and A^ is the adjoint operator of Nl^. If the operator is a potential operator, the potential 
F is given by 

F = JuJ^ N{Xu)dXdV. 

where / dV represents integration over the physical domain and dX represents inte- 
gration over the scalar variable A. We test the operator 

N{ip) = 2ipxxxt + Q^xx^xxxx + 69? 

XXX ~^ 'Pxxxxxx ~l~ (^'Pxxyy 

and find = N'^. Hence, 

Hf) = IvI^^N{Xip)dXdV 

= I V'lo [^^^xxxt + G^'^^xxfxxxx + (^^'^flxx + ^^xxxxxx + O- X(pxxyy]d^dV 
= J ^pi^Pxxxt + ^iflxx + 2ip 

XX '^xxxx "I" 2 'Pxxxxxx "I" 2 V^^yv 

To obtain F{u) in a more familiar form, integrate by part and discard the integration 
over the boundary since it has nothing with the Lagrangian density, then we get potential 

Fi'P) = J ifxx^xt - ^^Ixx + ^'Ply + flx)dV. 

So we can determined that the Lagrangian density for equation (2.2) is 

Hfi'P)) = ifxxfxt - \vlxx + '^'fly + 'Plx)dx AdyA dt. 
Corresponding to Lagrangian density L{j^{ip)), the action function is defined as following: 

'S'(¥') = / Lij^i'P)): J7 is a open set of X. 
Ju 

Let G be the Lie group of ttxy bundle automorphisms rfy covering r)x- Denote r)Y 
an smooth path in G such that 

Tp = r)^ o ip o {rj^y^ . 
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The vector field of 77y is 



yx 



We say that </? is a extremum of S is 

d , 



dX 



-.0 sm = 0. 



Now we consider the variation 



d_ 
dX 



a 



where 



X 




X 


V 




y 


t 




t 



A direct computation shows 



— |a=o S{ip) =h + h, 



in which 



XXX ^xxxxxx ~l~ (^fxxyy) 

{V^ - - ipyVy - iptV*)dx AdyA dt, 



(2.3) 



= lauK'PxxVxt - ^y^lxx + ^fly + flxW^ 

+iPxxiV;f - ^xV^ - ^yVy - iftV^ - ^xxV^ - ^xyVy - VxtV')]dx A dy 

[-{<PxxVxt - \^lxx + i'Ply + 'Plx)y^ 

+a^xy{V^ - ^xV^ - VyVy - VtV^ - ifixxV^ - ipxyVy - <PxtV')]dx A dt 

+ [{Vxx^xt - yixx + ^Vly + Vlx)y'' 

+ {V^t + SV^L + ^xxxx)iy^ - ^xV^ - fyVi - ^tV'x - VxxV- 
-fxyVy - iPxtV) - ^xxx{y$x - fxV^x - "^fxxV^ - yPyVi, 

-^'PxyVy - iptV^^ - 2ip^tV^ - ifxxyVy - (pxxtV^ - fxxxV') 

+ (-2(/Cxxt - ^^xx^xxx - (^^xyy - ^xxxxx){V'^ - ^x^ - ^yV^ - iptV*)]dy A dt. 

(2.4) 

In [5], the author introduced the Lepagean equivalents which is a generalization of 
the Poincare-Cartan form and proposed the Carton form is not necessarily unique for 

higher order. The author also pointed out that every Lagrangian density has a Lepagean 
equivalent on J'^^^'^Y . It turns out that Lepagean equivalents exists on jet bundles of 
order 2A; — 1 or higher, but not necessarily on jet bundle of lower order. The author 
presented that: As being the principal part of a Lepagean equivalents, the Cartan form 



5 



always exist and arc typically nonunique unless k = 1. we find it is caused by the higher- 
order mixed multiple integral in using Stokes' formula in actually calculus of variations. 
By I2, we can define a Cartan form 

X 'Pxxxxx 

-^xx^xt - '^vlx - Vxx^xxxx)dx Ady Adt 

+ {-2ipxxt - ^VxxVxxx — cr^xyy " ^xxxxx)d^p A dy A dt (2.5) 
+ {Vxt + ^^Ix + 'i^xxxx)dipx Ady Adt - ifxxxd^Pxx Ady Adt 
—a(pxydipx Adx Adt + ^Pxxd^Px Adx A dy. 



Since 



j^{tp)*dx = dx, j^{(p)*dy = dy, j^{(p)*dt = dt 
f (if)* dip = (fxdx + ifydy + ftdt, 
f{ip)*d(px = ^xxdx + ipxydy + ipxtdt, 
f{ip)*d(pxx = ^xxxdx + ipxxydy + ^xxtdt, 



we have 



JdU 



here j^iV) is the jet prolongation of the vector field V [8]. The multisymplectic form 
is the 4-form Jl^, = —dOz- Form 6l defines a multisymplectic structure on jet bundle 

J5(y). 

Now, we consider the Euler-Lagrange equation for the action function S{ip). 
Since L{j^{Jp)) =j^{ip)*9i^, we have 



- dX U=o J„A 



in which symbol £ denotes the Lie derivative. 
By the Cartan's Magic formula [6] 



We can obtain 
d 



-_oS{ip) = - [ f{ipnf{v)\nL)+ I j\^r{f{y)\eL) 

JU JdU 

If y is a vector field with compact support, we have 

JdU 

Hence, a necessary condition for ip to be an extremum is that 

JU 
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for any V with compact support. By compute the integral and obtain that 



+(TiPxxyy) {V^ - ^xV^ - ^y^^ - ftV) ^^'^^ 

Remark: We can get the Euler-lagrange equation from the vertical variation , and 
the and directions horizontal variations gives the law of conservation of momen- 
tum. The law of conservation of energy can be obtained along time-direction horizontal 
variation. 

Taking the TTxy-vertical vector field V and using the standard method from the 
calculus of variations, we obtain that satisfies 

'^^xxxt + ^^xx'^xxxx + 

XXX ^xxxxxx "I" ^^xxyy (^•'^) 

i.e., the equation (2.2). So, for any vector field V , 

f{^r{j\V)\^L) = G (2.8) 
holds. A short computation verifies that 

j\ipr{P\nL) = (2.9) 

where P € TJ^iY) and is T^^^g^^^ -vertical. For any W G TJ^{Y), there exists vector 
field V, such that 

W = f{V)+P (2.10) 

So, by (2.8)-(2.10), if if is an extremum of S, j^{ip)*{W\^L) mush vanish for any vector 
field W G rj^(y), thus, we get the Euler-lagrange equation 

/(<p)*(TyjOi) = (2.11) 

for any vector field W G TJ^{Y). In the following part, we consider the multisymplectic 
form formula and a corollary of the multisymplectic form formula. About the multisym- 
plectic form formula for first order field theories, please refer to [11]. 
Theorem 2.3 let Vjy and demote two one-parameter symmetry groups of equation 
(2.11) and the corresponding vector fields are V and W . Then we have the multisym- 
plectic form formula 

" f{vr{f{v)\f{wmL) = Q. (2.12) 



idU 

Proof: Since = [f{W),f{V)], from (2.11) we have 

= jujH'Pn[jHw),jHvmL) (2.13) 

Because rjy and .^y- are two one-parameter symmetry groups of equation (2.11), so for 
any vector field Q G TJ^{Y), we have 



^ \x=ojHri^o'Po{ri^xr'nQ\^L)=jHv^r£jHv){Q\^L) = 0, (2.14) 
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Thus (2.13) becomes 



x=ojH^YOVo{^^)-^y(Q\nL)=f{ipr£rHw){Q\^L) = 0. (2.15) 



= -jujH^nfiv)\£j^w)^L) 
= -luJ^y^nfivrnfmifiL)) (2.16) 

j^{V)\j^{W)\i^L can be written as 

f{v)\f{w)\nL = f{v)\difiw)\eL)-f{v)\£j^iw)0L 

= £j^v)(j'{W)\eL)-jHV)\£j^w)0L-dij^V)\j^W)\9L). 

(2.17) 

So from stokes' formula we can obtain that 

= IdujH^n£jHv)ifiw)\eL)-f{v)\£^s(^^)eL-d{f{v)\f{w)\eL)) 

= IujH^rd{£j5^v){fm\9L)-jHv)\£j^w)0L) 

= Jufi<pn£j5^y)£j5^w)0L + £j5^y)if{W)\nL)-d{fiV)\£j5^W)dL)). 

(2.18) 

By 

£jS^v)£jHw)OL = f{V)\d{£js^w)0L) + d{j5{V)\£j5^w)0L), (2.19) 

we have 

= Juf{v^njHv)\di£^5^w)0L) + £j5^v)ijHw)\r2L)). 

Hence, by (2.14) and (2.16), we obtain 



(2.20) 



dU 

Although the covariant Legendrc transformation ( or complete Legcndrc transforma- 
tion ) which transform time and space variable simultaneously are not necessarily unique, 
for this fixed Cartan form ^l, we can construct corresponding covariant Lengendre trans- 
formation of Lagrangian density L. Let v = <px,u = <Pxx,w = fxy^P = ^xt, taking the 
covariant Legendre transform of Lagrangian density L: 

~ "^Vxxt ^^xx^xxx (^Vxyy ^xxxxxi 
P^^ = Vxt + 3(/?L + "Pxxxx, P^* = Vxx, 

= aipxy, P^^^ = -(fxXX, 

According to the covariant Dc Bonder- Wcyl Hamilton function theories [18] and the 
multisymplectic concept introduced by Bridges [3], KP equation can be reformulated as 
a system of ten first-order partial differential equations which can be written in the form: 



MZt + Kl^ + Lly = VS'(Z), 



(2.21) 
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where 



M = 



( ° 
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S{Z) =up + ^ {p^^^f + ^w'^ + - p^'v - p^'^'u - p^*p - p^^^w. 

V/S is the gradient of S with respect to the standard inner product on M.^'^. The 
system (2.21) is a Hamiltonian formulation of the KP equation on a multisymplectic 
structure. Although this formulation is not the best, we give a constructive method to 
get multisymplectic form which proposed by Bridges [3]. It should be pay attention that 
(2.21) is different with the multisymplectic form of KP equation in [4]. For equation 
(2.21), there is a conservation law 

dt{dZ A MdZ) + da:{dZ A KdZ) + dy{dZ A LdZ) = 0. (2.22) 
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Substituting M,K,L into (2.22) leads to 



§l{d(px A difxx) + -^{dy^ A {-Gcpxxxdifxx - ^Vxxd^Pxxx — cfdipxyy — 2d(pxxt - difxxxxx) 
+dipx A {Qipxxdifxx + dipxt + dipxxxx) + dipxx A d{-ipxxx)) + i^id^x A d{aipxy)) = 0. 

(2.23) 

This multisymplectic conservation law (2.22) consistent with our theorem 2.3. We 
can regard the conservation law (2.23) as a corollary of the theorem 2.3. Let V,W he 
TTxy-vertical and have the expressions V'^-^,W'^^. Thus the corresponding j^(y) and 
fiW) have the coordinate expressions {V^, , F/, , V^,, V^%, V,t, V^y, V^^, V;f,,, 

tyy^ 



yV Y<fi Y<p Y<fi V'^ V'^ V'^ 

xtx ' xxy ' ' xtt ' ' xyy i ' ttt' yyy ' tyy ' tty ' xxxx ' ' xxyy ' * xxtt ' xxxy ) ' xxxt i ' xvvv ' xttt i ^ tttt ' 
Y<p Yf yf Y<p y<p yP y^ yf yp yf 

' til It II I "I ' Tt It 111 t » TTI I'll 1 ' TTTt It ''Y"Y">"Y»'Y»1 ' -y-yn^n^tl t " 'Y"Y"Y"Y»T 1 " -T"! CI I 1 " 'yy' i^TT •> ' •Wtltltl', »'> 



xxxt 1 ' xyyy 
' xxttt ' xyyyy 



xttt' ''tttt' 
Vf„„„„, Vxtttt' 



xtx ' '' xxy ' xtt ' '' xyy ' " ttt' " yyy ' '' tyy ' '' tty ' 

vy" vjf * 

' ' tyyy ' ' ttyy ' ' ttty ' xxxxx ' xxxxy ' xxxxt ' xxxyy ) 

^yyyyyi ^ttttti ^yyyyti ^yyytt' ^yyttt' ^ytttt) ^^'^ (VF''', W!^, , , W^x' ^yy> ^tt ' ^xt' ^xy' 

ly"^ VF'^ T^'^,, W^t Wf Wft W'^ 

" " xxx' " xtxi "xxy " xtti "xyy "tttt "yyy "tyy "tty " xxxxi "xxyy "xxttt "xxxy 



ytt 



yyy 



tyy "tty "xxxxt 



, Wf,,, Wf Wf wf, W^ Wf W^ , Wf 

" xxxtt " xyyy t " xtttt " ttttt " yyyy t " tyyyt " ttyyt "tttyt " xxxxx t " xxxxy t " xxxxt t " xxxyy t 



w. 



,, Wf 

xxxtt t ' ' xxyyy 



W"^ ,„ Wf ' ' ' W"^,,,, Wf 

' ' xxttt t ' ' xvvvv t ' ' xtttt t " vu 



xyyyy t 



yyyyy 



Wf„„ W^ , W 

tttttt yyyytt "^"^i 



xxxxy t 

f wf wf 

yyyttt ^^yytttt ^^ytttt)- 



We can compute 



jH^njHv)\jHw)\nL) = (wf^vf - wfvf)dx a dy 

+ - '^WfyVf)dx Adt + [Wf,Vf, - Wf^^Vf, 

-wf{v:^t + ^v^^vf, + yfxxx) + y'f(wtt + ^^xxWf^ + w^,, 

+Vf{-2W^,t - 6^xxWf,, - 6^,,,Wf, - aWf - Wf,,,,) 



+W^{2V:^^^ + Q^, 



+ Vf 

xyy ' ^ xxxxx 



)]dy A dt. 



By the Stokes' formula, we can obtain 



Since U is arbitrary, we get 

I - W;fVf ) - i^iaWfVf - aWfyV^) + 9_[Wf,Vif„- 

W^xV^ - WnV^t + ^^xxV^x + Vitxxx) + VxHW^t + ^VxxW^x + W^xxx) 

^Vf{-2Wt,t - Q^xxWf,, - Q^xxxWf, - aWfyy - W^,,,,) 

+W^{2Vl, + Q^xxV^xx + Q^xxxV^x + ^V^yy + V^xxxx)] = 0- 



(2.24) 



If we let 
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M\x, y) = x'^M^y, W^x, y) = x^M^y, M^ix, y) = x'^M^y (2.25) 

and set 

n^(v) = (v^ v'^, ) 

J y I V' ■> ^ X 1 " XXI ' xyi " xti " XXX i ' xyy i " xxt ' ' xxxx i ' xxxxx ) 

the others coordinates vanished. Then conservation law (2.24) can be written as 

^^M\f{V),f{W)) + ^M-if{V),fiW)) + d_My(j5(v),jHW)) = (2.26) 

Also since the translation invariance of the KP equation, we choose V = W = and 
take it into (2.26), conservation law (2.26) becomes: 

§i{dipx A difxx) + -§i{d^ /\ {-^(Pxxxd(Pxx - QfxxdiPxxx - Crdfxyy - 'idipxxt - dfxjcxxx) 

+dipx A {Gcpxxdifxx + dipxt + dipxxxx) + dipxx A d{-ipxxx)) + §^{d^Px A d{aipxy)) = 0. 

(2.27) 

i.e. the conservation law (2.23). 

In the numerical study, the multisymplectic conservation law can be used to design 
multisymplectic numerical schemes, i.e., numerical schemes which can preserve the mul- 
tisymplectic conservation law. 
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3 Multisymplectic Preissman Scheme for the KP Equation 

In this section, we consider the multisymplectic Preissman scheme for the KP equation. 
The equation (2.21) can be reformulated as 



dx ~ ^' 



dx ' "Sy ' dt 

aw — p^y = 0, 

„xt 



dp' 



+ 



dp^ 



-p 



U — p 

dip 



0, 



(3.1) 



w, 



ax 
dv 
dy 

dt 

_ ^ pXXX 



For convenience, we assumed that the spacing of the grid points in the x, y, t directions 
is uniform respectively. We apply the implicit midpoint discretization in time and in space 
to (3.1), and obtain that: 



Ax 

Ax 

xt 



0, 



+ 



xy xy 
i+l,j+l,k+l i+l,j,k+l 



Ay 



At 



i+l,j+^,k+^ i,j+^,k+^ 



pXX ^ ^ ^ 

»+2 iJ+2 ''''+2 



"^^1+ 2 J+ 2 1^^+ 2 



'^i+l,j+^,k+^~ 



(3.2) 



Ax 



^i+l,j+^,fc+^" 



'"i+^,j+l,k+^ 



Ax 
1 



Ay 



w, 



1 , 



V T ^ 



At 

^■11 ■ I 1 I I 1 ~" 



Ax 



Where Ax is the x-direction step. Ay is the y-direction step. At is the time step and 
""i+i j+i fe+l = u{iAx + ^,jAy + ^,kAt + ^), the others is similar. 
In fact, the discretization result lead to the pressman scheme 



i,M(Z, 



At 



'2 



(3.3) 
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The (3.3) preserve the discrete multisymplectic conservation law: 



dv 



-dv 



Adp^ 



-+ 



»+7,J,fe+' 



_ Adp 



Adp^ 



Ay 



+I,j+i2,k+i2 ^ij+^,k+^ "'-ij+i,fc+^ 

Ax 



(3.4) 



du. 



Ax 



1 — du. . , 1 , , 1 Adp- , , 



Ax 







Although the Preissman scheme (3.3) is multisymplectic, it took too efforts to realize. 

Hence we elimate the auxiliary variables ip,v,w,p,p^ ,p^^ ,p^y ,p^*,p^^^ by a trivial com- 
putation and obtain the following multisymplectic forty-five points scheme: 



2^2( k 



-I- 6n: 



} 



)} 

where we denote ^ = u{iAx,jAy,kAt), f^j = 3{u{iAx,jAy,kAt))'^ and 



(3.5) 



« J ' y « J 
,fc I 



^A; _i_ fk 

■'i+hj i- 



' 2 'J 



(3.6) 



4 Some Numerical Results on Soliton and Solitary Waves 

In this section, we test the forty-five points scheme on soliton and solitary waves over 
long time intervals. 



0- ...^^ 




i.. 



™' 3/ 7 



7 'T° 



Fig.l: one soliton at t=0. 



Fig.2: one soliton at t=10. 
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At first, one line soliton of KPI equation is considered. We choose small interval on 
y-direction just for computing convenience, it have nothing with the scheme and result. 
We take the following initial conditions 



a/2 

u{x,y,0) = 2sech^ {x —y — 6) 

and the exact boundary condition. The KPI equation has the theoretic solution 



(4.1) 



u{x, y, t) = 2sech (x 



V2 



-y 



6) 



which represents one line-soliton propagating with the velocity | in the direction with the 
angle of tan~^{\/^) to the positive x-axis. We take the test on the domain [0, 40] x [0, 2] 
and choose Ax = 0.2, Ay = 0.1, At = 0.01. Figure 1 shows the initial condition and 
Figure 2 shows the numerical solution at time t = 10 . In fact, it propagating 25-unit 
distant as we indicated. 

Next wc try two linc-soliton interaction. 

We take the initial condition 



u{x, y,0) = 2 ^ kfsech?[ki{x + A^y - xo,i)] 
=1 



(4.2) 



where k\ = 1.0,^2 = l/v^jAi = — l/-\/3,A2 = —1.0 and Xq.i = 6.0,xo.2 = H-O. and the 
exact boundary condition. The initial condition (4.2) corresponds to two line-solitons, 
each with amplitude 2/c,? placed initially at x = xo,i and moving with velocity vi = 
Akf — 3A? along the x-axis (z=l,2). 





15 l/ 



Fig.3: soliton interact at t=0. 



Fig.4: soliton interact at t=1.5. 



We carried out the computation on the domain [0,40] x [0,2], and choose Ax = 
0.2, Ay = 0.1, At = 0.01. The initial condition (4.2) is showed in Figure 3. The larger 
line-soliton on the right will moves with a velocity 3.0 to the positive x-direction and the 
smaller one on the left will moves with a velocity 1.0 to inverse direction. As time go 
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on, they will collide with each other, it is show in Figure 4. Figure 5 shows that the two 
line-solitons have separated completely after collide and restored their original shape by 
the t = 3. 

Fig.5: soliton interact at t=3. 
We also consider the lump type solitary waves of the KPI equation. The lump type 
initial condition used for test is 



(-{x-xof + fiHy-yo? + ji 

n(x,y,0)=4- 



((x - xq^ + n'^{y - + 



2 I 1^ ' 



(4.3) 



where the parameters )U^=1.0, xo=10.0, yo=10.0, and we adopt the exact boundary 
condition. 



. . .....^ ^, . ... 



.-1 i<!i!!i iiJE. ■■ ^T^'?^■s^~.,iv.■.■.;^:■■■;^^^ 



" " IB 14 12 ,/ 
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Fig.6: lump type solitary wave at t=0, Fig.7: lump type solitary wave at t=0.5. 

We compute in a rectangle [0, 20] x [0, 20] and choose Ax = 0.1,A?/ = 0.2,At = 0.01. 

Figure 6 shows the initial condition and figure 7 shows the numerical solution at time 
t = 0.5. By the time t = 1, the result is showed in figure 8. The lump solution of the 
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KPI equation can be expressed as 

{-{x ~xo- SfiHy + i^^iy - yof + ^ ) 



u{x,y,t) = 4- 



(4.4) 



According to (4.4), this lump type solitary wave will move to the positive x-direction 
with velocity 3//^. We can see the moving of the lump solitary wave from the graph. 



*^'£-!SK - Li-.^'*r^ ■ 
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Fig.8: lump type solitary wave at t=l. 

A future task is expected to find a proper numerical boundary conditions for collision 
of the two lump type solitary waves. 
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